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Abstract 

We provide a class of QFTs which exhibit dissipation above a threshold energy, thereby 
breaking Lorentz invariance. Unitarity is preserved by coupling the fields to additional 
degrees of freedom (heavy fields) which introduce the rest frame. Using the Equivalence 
Principle, we define these theories in arbitrary curved spacetime. We then confront the 
trans-Planckian question of inflationary cosmology. When dissipation increases with the 
energy, the quantum field describing adiabatic perturbations is completely damped at the 
onset of inflation. However it still exists as a composite operator made with the additional 
fields. And when these are in their ground state, the standard power spectrum obtains 
if the threshold energy is much larger that the Hubble parameter. In fact, as the energy 
redshifts below the threshold, the composite operator behaves as if it were a free field 
endowed with standard vacuum fluctuations. The relationship between our models and 
the Brane World scenarios studied by Libanov and Rubakov displaying similar effects is 
discussed. The signatures of dissipation will be studied in a forthcoming paper. 
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1 Introduction and presentation of the settings 



Even though relativistic QFT provides an excellent description of particle physics, being non- 
compact, Lorentz symmetry can only be tested up to a certain energy scale. Thus one cannot 
exclude that some unknown high energy processes break the invariance under boosts, thereby 
introducing a threshold energy Alv, and a preferred frame. It is therefore of interest to deter- 
mine what would be the signatures when this possibility is realized pp. 

This question is particularly relevant for inflationary cosmology because primordial density 
fluctuations arise from vacuum fluctuations which had very short wave lengths (very large 
proper frequencies) at the onset of inflation[2]. Indeed, the initial frequencies of the modes we 
observe today in the CMB anisotropies were all larger than 

n in = He N ^% (1) 

where H is the value of the Hubble parameter during inflation, and where N extra is the number 
of e-folds from the onset of inflation till to, the moment when the (comoving) scale of our 
visible universe exited the Hubble radius. (The total number of e-folds is thus N + N extra 
where N ~ 60 is the number of e-folds from t till the end of inflation.) Irrespectively of the 
inflationary scenario, VL in is thus larger than the Planck mass Mpi when N extra > In Mpi/H. 
Since H should be of the order of 10~ 5 Mp;, the initial frequencies of observable modes were all 
trans-Planckian when N extra > 5 In 10 ~ 12. Moreover, since in most scenarios Nextra ^ 12, 
Qi n was generically much larger than the Planck mass. 

In the absence of Quantum Gravity, there is no understanding of the nature of the degrees of 
freedom at these scales: we neither know what is their Hilbert space, nor how they propagate, 
and this, because the notions of differential geometry used when describing (quantum) fields 
in expanding universes might not apply at very short distances. In particular eq. (CQ) which is 
based on a relativistic mode equation in an accelerating universe might loose its validity above 
a certain Ultra- Violet scale Alv- Since this scale might significantly differ from the Planck 
mass, we shall keep it distinct from the Planck mass. We shall also use the abbreviations LI 
and LV for Lorentz Invariance and Lorentz Violation respectively. 

In this context, it is instructive to phenomenologically parameterize the deviations from LI 
and determine their signatures on inflationary spectra. This line of thought has been proposed 
in the context of Hawking radiation from black holes, wherein the asymptotic quanta also arise 
from vacuum fluctuations with exponentially growing frequencies [3111] . The deviations from LI 
have been characterized by non-linear dispersion relations imposed in a particular rest frame, 

Q 2 = F n {p 2 ) = p 2 (1 ± p 2n /A% + 0(p 2n+2 )), (2) 

where Q and p are respectively the frequency and the norm of the spatial momentum measured 
in that frame. The first deviation is characterized by a power of p/ Alv and the sign defining 
super-luminous (+) and sub-luminous (— ) cases. When Alv is much larger than Hawking 
temperature, and when the frame is freely falling, it was shown that the asymptotic properties 
of Hawking radiation are unmodified, even though the near horizon propagation is radically 
modified when Q is larger than Alv- The robustness of asymptotic properties relies on the 
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adiabatic character of the evolution of the vacuum (ground) state [5]. The same program was 
then applied to the inflationary spectra [HllTj, and similar results were obtained because, in this 
case as well, the ground state adiabatically evolves when a = H/Alv <C 1 [8]. However, so far, 
only dispersive models have been systematically studied. 

The aim of the present work is to extend this analysis to dissipative models. To this end, 
we shall first construct QFT displaying dissipative effects in the UV. Indeed, unlike dispersion, 
dissipation requires enlarged dynamical settings because if one tries to introduce dissipation 
from the outset in eq. (T2]), one looses both unitarity and predictability. To preserve them, we 
shall work with Hamiltonian theories in which dissipative effects are caused by interactions with 
additional fields. Doing so, we shall establish that dissipative effects are generic. That is, when 
starting with a Lagrangian in which LI is broken in the UV, the effective theory unavoidably 
develops dissipation above a certain energy scale, simply because nothing can prevent this. 
(With relativistic QFT instead, LI did prevent it). 

Since we want to construct generalized QFT, we should decide what to keep. First, we want 
to preserve the unitarity of the description because the calculation of the power spectrum of 
(adiabatic) density fluctuations [2] rests on the identification of a scalar field, hereafter called 
0, which obeys Equal Time Commutators (ETC). This identification is necessary to fix the rms 
amplitude of density fluctuations in the vacuum. Explicitely the power spectrum is given 

P p (t) = 4V [ (P) 3 e^G a (t,x;t,0), (3) 
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where G a is the anti-commutator of evaluated in the asymptotic (Bunch-Davies) vacuum 
much after horizon exit. In the absence of dissipation, P p can be deduced from the norm of 
<f) p n (t), the Fourier mode of with asymptotic positive frequency. However in the presence 
of dissipation, this notion of (free) mode disappears. Hence the knowledge of G a becomes 
necessary since eq. (J3]) gives the only way to obtain the power spectrum. 

To get dissipation we shall thus introduce additional degrees of freedom, hereafter called 
Then the whole system (0 + \]/) will evolve unitarily, by construction. This guarantees that the 
ETC of is exactly preserved (in a non trivial way since undergoes dissipation). Moreover, 
the (dressed) 2pt function of will be given by the usual QM trace 

G w (x,y) = AR[pT$(x)4>(y)], (4) 

where is the initial matrix density of the total system, where 0(x) is the Heisenberg field 
operator evolved with the time ordered exponential of the total Hamiltonian, and where the 
trace is taken over both \l/ and 0. The anti-commutator G a determining the power spectrum 
in eq. (jSJ) is simply the symmetric part of Gw- 

Our second requirement concerns the properties of dissipative effects. When considering 
the theory in vacuum and in Minkowski space-time, we impose that these effects preserve 
the stationary, the homogeneity, and isotropy of flat space-time. These requirements define 
a preferred frame which is inertial and globally defined. Then irrespectively of the choice of 
\l/ and \l/-0 interactions, the Fourier transform of the retarded Green function, G r (x,y) = 
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&(t x — ty) 2lmGw{x, y), is of the form 



G r (uj,p) = r- . (5) 

(-^ +J 9 2 + £ r (^)J 

In the true vacuum, at the level of the 2 pt functions, the dissipative (dispersive) effects are 
thus completely characterized by the imaginary and odd part in uj (real and even part) of the 
self energy Ti r (uj,p). 

In these expressions, the energy u and the momentum square p 2 are defined in the preferred 
frame. The novelty is that S r is a function of uj and p separately, and not only a function of 
the relativistic invariant uj 2 — p 2 as it is the case in LI QFT. Therefore dissipation can become 
significant above a critical energy on the mass shell, i.e. along the minima of the denominator 
of eq. (J5]); a possibility forbidden in LI theories. For instance, one verifies that the following 
self-energies induce significant dissipation only above A^y 

(2 \ n 

With this equation we have identified the relevant (lowest order) quantity governing dissipation, 
the equivalent of the first order deviation in eq. In the body of the paper, we shall provide 
Lagrangians of \l/ and which give this class of self-energies labeled by n. Rather than focusing 
on a particular case, we shall describe the whole class of dissipative behaviors (at the level 
of 2-point functions). We shall thus follow a phenomenological approach based on minimal 
assumptions (i.e. unitarity), rather than some "inspired" approach (e.g. by string theory, or 
branes scenarios P, EH]) which would lead to a particular sub-class of models. 

Our third requirement concerns the extension of our models from Minkowski space to curved 
space-times. To define our QFT in curved space, we simply implement the Equivalence Principle 
(EP). It fixes the action density our models to be the covariantized version of that we had in 
Minkowski space (up to some non-minimal coupling). To perform the covariantization, it is 
useful to characterize the preferred frame in a coordinate invariant way by P, a unit time-like 
vector field [TT] . In terms of this vector, u and p 2 are given by 

u = Fp^ p 2 =^ PllPu , (7) 

where _!_'"'= + l^V is the (positive definite) metric in the spatial sections orthogonal to Z M . 
The covariant action will be a sum of scalar functions of the four local fields 0, g^ v and l v 
which reduces to the Minkowski model in the zero curvature limit. Very importantly for us, we 
shall see that the EP guarantees the adiabaticity of the evolution of the (interacting) ground 
state as long the gradients of the metric are much smaller than the UV scale A^y. 

For the interested reader, we mention that additional comments which place the present 
work in broader contexts can be found in the introduction of [TJ]. These comments concern 
Quantum Gravity and the description of " mode creation" in expanding universes. 



4 



2 Dissipation in Minkowski space from LV effects 



In this Section, we provide a class of models defined in Minkowski spacetime which exhibit 
dissipative effects above a certain energy scale A LV . Stationarity, homogeneity and isotropy 
will be exactly preserved. Therefore, the only invariance of relativist QFT which is broken is 
that under boosts. These theories define a preferred rest frame which is globally defined, as it 
is the case in FLRW space-times. For the simplicity of the presentation, we shall first work in 
the preferred frame. At the end of this Section we shall covariantize our action and generalize 
it to curved spacetime. 

2.1 Free field settings 

We start with a brief presentation of the free field quantization to introduce notations and to 
point out what are the properties which will be lost in the presence of dissipation. 
The action of our free massless field is the usual one: 

S<t> = \f dtd 3 x (d t <p 2 - «9 X • «9 X 0) , (8) 

where t and x are Cartesian coordinates. Due to the homogeneity of space, the equation of 
motion can be analyzed mode by mode: 

^' x) = /(2^ eiPX0p(t) - (9) 

The Fourier mode 4> p (t) obeys 

(^ + ^ 2 )0 P = O, (10) 

where u 2 = p 2 = p • p is the standard relativistic dispersion relation. Notice that eq. ( fTOl) is 
second order, homogeneous (no source term), and time reversible (no odd power of d t ), three 
properties we shall loose when introducing interactions breaking LI. 

In homogeneous space-times, the Equal Time Commutator between the Heisenberg field 
operator and its momentum implies that the mode operator (q-number) 4> p (t) obeys 

When decomposing this operator as 

p (t) = a p p (t) + a t _ p 0;(t), (12) 
where the destruction and creation operators satisfy the usual commutators 

[«P> a p'] = <5 3 (P-P')> [a p ,a p /]=0, (13) 
eq. fllip is verified because the Wronskian of the positive frequency (c-number) mode 

P (t)=e"^7(2^) 1/2 , (14) 
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is constant (and conventionally taken to be unity). 

Had an odd term like 7^ be present in eq. (II 01) the constancy of the Wronskian would have 
been lost. Hence the possibility of realizing the ETC (TTTj) with the help of eq. ( 1131) would have 
been lost as well. This already indicates that, unlike dispersive (real) effects, dissipative effects 
require more general settings than the above. 

2.2 Interacting models breaking LI, general properties 

We now introduce additional degrees of freedom, here after collectively named which induce 
dissipation above the energy A^y. We shall work with a particular class of models in order to get 
an exact (non-perturbative) expression for the two-point function of eq. (jl]). Before introducing 
these models, we derive general results valid for all unitary QFT's possessing dissipative effects 
above A^y in the ground state (the interacting vacuum). 
We assume that the total action decomposes as 

St — S<p + S-ip + Scf,^ , (15) 

where the first action is that of eq. ([8]), the second one governs the evolution of the fields, 
and the last one the coupling between and these new fields. We also impose that the last two 
actions preserve the homogeneity and isotropy of Minkowski space but break the invariance 
under boosts. From now on, the coordinates t, x are at rest with respect to the preferred frame 
defined by + S^,*, i.e., d t = l^d^ in the covariant notation of eq. (J7j). 

When the state of such system is homogeneous, the Fourier transform of the 2pt function 
of eq. (jl]) is of the form 

G P;P ,(t,t') = Tr[p T j> p (t)4>l,(t')}, 

= G w (t,t';p)5 3 (p-p'). (16) 

At this point, an important remark should be made. In the presence of interactions, the notion 
(and the usefulness) of the time-dependent modes of eq. ( fill) has disappeared whereas the 
function Gw(t,t';p) of eq. (ITBl) is always well-defined, for all choices of Hf and + S^^. 

When the situation is stationary, Gw(p] t, t') further simplifies in the frequency representa- 
tion: 

Gw(t,t';p) = J ^e-^-^G w (u;,p) . (17) 

When working in the ground state, one finds the simplest case, because whatever the \I/ fields 
may be, the Fourier transform of the time ordered (Feynman) propagator, 

2G F = ReG w + i lmG w sign(t - tf), (18) 

is always of the same form as G r in eq. and therefore characterized by a single (complex) 
function T^p^Ujp). When restricting attention to Gaussian models, Y^pioj^p) is given by a 1-loop 
calculation, whereas in general, it contains a series of 1PI graphs. 
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In non-vacuum states and in non-stationary situations, the 2pt functions and self-energies 
have a more complicated structure [20]. To understand this structure it is useful to study 
separately the commutator G c = ilmGw and the anti-commutator G a = ReGw- 

Let us conclude with two remarks. First, the effective dispersion relation of <fi is a posteriori 
defined by the poles of eq. (j3J) which are governed by T,(u,p) [131 [T3]. In this way, non- 
trivial dispersion relations arise from dynamical processes rather than from being introduced 
from the outset. The present work therefore provides physical foundations (and restrictions, 
as later discussed) to the kinematical approach which is usually adopted pQ. Second, from 
analyzing dynamical models, we shall see that, even in the vacuum, on-shell dissipative effects 
(i.e. dissipation arising along the minima of the denominator of eq. (131)) are unavoidable when 
LI is broken in the UV by the action S& + S^, in complete opposition with the fact that 
on-shell dissipation is forbidden when working (in the vacuum) with LI actions. 



2.3 Gaussian models 

To simplify the calculation of Ti(ou,p) and to get non perturbative expressions, we assume that 
the action St is quadratic in all field variables. At first sight, this could be considered as an 
artificial hypothesis. However, it should be recalled that we are not after computing £ from 
first principles. Rather we aim to compute the signatures of the power spectrum (J31) given the 
properties of E, following the approach adopted in [31 El [7] . 

Given that we are preserving the homogeneity of Minkowski space, the Gaussian assumption 
implies that the total action splits as 

S T = Jd 3 pS T (p), (19) 

where each action Sr(p) depends only on P and the p-th Fourier component of The 
structure of these actions is 

s t (p) = \ f dt <t>;(-d? - + ^ J dt **( P )(-^ 2 - n?(p))*i(p) 

+£, j dt gi (p)<f) p y*(p), (20) 

where % is a discrete (or continuous) index, where Qi(p) is the energy of the quanta of the 
oscillators ^(p), and where gi(p) is the coupling constant at fixed p, i. Since P = 0_ p , 
^J(p) = 1 S r i(— p), St(p) + St(— p) is real. Instead of choosing a priori one model, it is more 
instructive to solve the equations of motion without specifying the set of ^j(p), their energy 
f£j(p) and their coupling gi(p). We shall choose them in due course. To preserve stationarity 
in Minkowski space, the fi's and the g's must be time independent. When &f{p) ^ Mf +p 2 , 
the kinetic action of ^(p) breaks LI and defines the preferred frame. On the contrary when 
^1 (p) = Ml + p 2 the preferred frame is only defined by S^ through the p-dependence of the 
coupling functions gi(p). 

Models of this type have been used for different purposes. They have been introduced 
(in their continuous version) to study non-pertubatively atomic transitions see [IS] and refs. 
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therein, see also [161 [17] for an application to the Unruh effect. They have been used in Quan- 
tum Optics [18], to model quantum Brownian motion, and to study decoherence effects [19J. 
Depending on the aims, they can be solved and analyzed by means of different methods. In 
what follows we shall use the simplest approach based on Heisenberg picture^] Since we shall 
work at zero temperature, this approach offers a simple characterization of the state of the 
system in terms of the ground states of the free modes before <f> — ^ interactions are turned on. 

To prepare the application to inflation, we shall treat u 2 , Vtj and gi as arbitrary functions 
of time (in cosmology these quantities become time dependent through their dependence in the 
scale factor a(t)). The equations of motion are 

(^ 2 + ^)0 p = S^ i (p)vI/ J (p), (21) 
(d? + tf)y i (p)=g l (p)<p p . (22) 

The general solution of the second equation reads 

* i (p,t) = V?(p,t) + J dt f R°(t,t';p)g l (t';p) ( J )p (t'), (23) 

where $?°(p,t) is a free solution which depends on initial conditions imposed on \I/i(p). The 
second term contains R°(t,t';p), the (free) retarded Green function of (p) . It obeys 

(d? + tf(p))R°(t,t>;p) = 5(t-t>), (24) 

and vanishes for t < t'. Injecting eq. (J21I1) in eq. f[2"TT) one gets 

(^ 2 + ^)0 p = S l ^(t;p)vl>°(p,t) + S i ^(t;p) J dt'R^t'-pMt'-p)^'). (25) 

The solution of this equation has always the following structure 

= + / dt'G r {t,t'- 1 p)[T li g i {t' ] p)^° i {p 1 t')] . (26) 

The first term is the "decaying" solution. It contains all the information about the initial 
condition of p . The second term is the "driven" solution. It is governed by the initial conditions 
of ^j(p) and by the (dressed) retarded Green function, the solution of 



/ 



dt'[5(t - t') (<9 2 + u 2 p ) - 9l (t; p) R°(t, t'; p) 9i (t'; p)} G r {t', h; p) = 6(t - h) . (27) 

Notice that 0j(t') is an homogeneous solution of this equation. Therefore the evolution of 
both (f) d and G r fully takes into account, through the non-local term in the above bracket, the 



1 Even though legitimate, we shall not use the general methods (Influence Functional, Master Equation, 
Closed Time Path Integral) [T5] which have been developed to study "open quantum systems" because they 
somehow hide the simplicity of the present models. Moreover we are planning (in a subsequent work) to study 
the correlations between (j> and . Therefore we shall treat <j> and on equal footing as in [15l [TBI HZ] • 
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back-reaction due to the coupling to the In Gaussian models, it is quadratic in ^. Hence 
(ft d and G r are series containing all powers of gi. Moreover since g%{t\p) are arbitrary functions 
of p and t, at this point, there is no reason to consider non-Gaussian models. 

To conclude this subsection, we notice that eq. (12"6"1) also furnishes the exact solution for the 
(Heisenberg) mode operator <f> p (t) because the equations we solved were all linear. Since the 
power spectrum in inflation is obtained from vacuum fluctuations, instead of further analyzing 
the time dependence of the mode (as one would do in classical terms), it is more relevant to 
study the correlation functions of (ft p (t). 

2.4 Structure of two-point correlation functions 

This sub-section mainly contains well-known results which follow from the linearity of eq. ( |26l) . 
The key result we shall later use is given in eq. (}3~Tj) . 

Since our models are Gaussian, the 2pt function of eq. (TT01) governs all observables built 
with the Heisenberg field operator (ft. To analyse it, as already mentioned, it is appropriate to 
study separately the commutator and the anti-commutator. 

We start with the simple part, the commutator 

G c (t,t';p)6 3 (p-pi) = Tr[p T [0 p (t),^,(t')]-]. (28) 

From eq. ( 1261) one sees that it decomposes into two terms, one due to the non-commuting 
character of <ft d , the other due to that of In addition, since both commutators are c- 
numbers, it is independent of px, the state of the system. Hence, for all Gaussian models, one 
has 

G c (t,t';p) = [<ft d (t),<ft d (t')}^ + J J dt x dt 2 G r {tM)Gr{t',t2)D{t u t 2 ) , (29) 

where the " dissipative" kernel D(ti,t 2 ) is given by 

D(h,t 2 ) = ^g^g^h) [*?(tx), *2(* 2 )]_ = ^g^G^^g^h). (30) 

Notice how this kernel combines the coupling g t and the non-commuting properties of 

The next property of G c is more relevant. To all orders in ^ and for all sets of g^ (even 
with arbitrary time dependence), one obtains 

id t G c (t,t';p)\ t=t , = l. (31) 

This identity corresponds to the ETC of eq. ( flTT) . The 1 on the rhs is guaranteed by the 
Hamiltonian character of the evolution of the entire system (ft + It is therefore this equation 
which replaces the constancy of the Wronskian that was relevant in the case of free evolution. 

Eq. (13~T|) is crucial for us for two reasons. First, since the operator (ft d {t) in eq. ( }2"6"1) 
decays in t — t in where t in is the moment when the interactions are turn on -because it is an 
homogeneous solution of eq. (1271)- the first term in eq. ( J29l) decays as exp — y(t + t' — 2t in ). 
Therefore when j(t + 1' — 2tj n ) ^> 1, the non-commuting properties <ft p are entirely due to those 
of the environment degrees of freedom, Secondly, these sum up exactly to 1, as if the driven 
term of eq. (121)]) were a canonical degree of freedom. 
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We now analyze the anti-commutator, 



G a (t, H- p) 5 3 (p - p') = Tr[p T {0 p (t), ^(f)}+ ] • 



(32) 



When the (initial) density matrix factorizes, = p<i>, as it is the case in the "free" vacuum 
before the interactions are turned on, G a also splits into two terms, 



The first term only depends on the initial state of (J). Similarly, the driven term only depends 
on the state of the environment through the " noise" kernel 



As for the commutator, in the presence of dissipation, the first term exponentially decays, 
expressing the progressive erasing on the information contained in the initial state of <fi. At 
late times therefore it is the state of \I/ which fixes the anti-commutator of <fi. This allows to 
remove the restriction that initially the density matrices factorizes. If one is interested by the 
late time behavior, only iV matters. 

In brief, we have recalled two important results. First, at late time, the Heisenberg field 
(f) reduces to its driven term, the second term of eq. fl2T|) . since both its commutator and 
anti-commutator are determined by those of $?°. Second, only two (real) kernels determined 
by the environment govern the two-point functions of 0, namely D and N of eqs. ( |30l l34"l) . 
Therefore the set of environments (Gaussian or not Gaussian) possessing the same kernels will 
give rise to the same 2pt functions. Hence they should be viewed as forming an equivalent 
class. The degeneracy can be lifted by considering correlations with observables containing the 
operators or higher order correlations functions of (for non-Gaussian environments), two 
possibilities we shall not discuss in this paper. 

To compute G c and G a , two different routes can be adopted. When gi(p) and fl(p) are 
constant, one should work in Fourier transform because the equations can be algebraically 
solved, in full generality. Instead when gi(p) and/or Q(p) are time-dependent, as it will be the 
case in expanding universes and in curved space-times, it is appropriate (but not necessary) to 
exploit the above mentioned degeneracy by choosing the and their frequencies fif (p) so as 
to simplify the time dependence of the equations. In the text, we proceed with time dependent 
approach. In Appendix A, we present the Fourier analysis which is straightforward. We invite 
the reader unfamiliar with the Quantum Mechanical treatment of dissipation to read it. 

2.5 Time dependent settings 

In Appendix B, we provided a class of models characterized by the power of p/Aiv which 
specifies how dissipative effects grow, see eq. (}85l) . This class covers the general case and can 
be used as a template to study the consequences of dissipative effects. In addition, as noticed 
after eq. ( |76l) . the dissipative effects are governed by D = HigfR°. Therefore all environments 
delivering the same kernel give rise to the same (stationary) phenomenology. 



G a (t,t';p) = Tr[p {0 d (t) ,0 d (t')}+] + // dt x tt 2 G r {t,t x )G r W \h)N{t u t 2 ) ■ 




(33) 



N(h,t 2 ) = E^Trfp* {gfaWfa) ,<fc(t 2 )^(t 2 )} + ] . 



(34) 
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In this Section we exploit this freedom having in mind the transposition of our model from 
Minkowski space to cosmological, and thus time-dependent, metrics. Therefore the selected 
models should possess two-point functions with simple properties when expressed in the terms 
of time (as opposed to Fourier components). The core of the problem is that G a depends, see 
eq. (IHBl . on the retarded Green function of <p which is not known. Indeed G r is only implicitly 
defined as a solution of eq. ( 1271) which is, in general, a non-local differential equation. We 
are thus led to choose \1/ in order for this equation to be local. This implies that the retarded 
Green function of \1/ appearing in eq. ( 1271) be proportional to S(t — t'). (Since this requirement 
concerns the Green function of the environment it does not restrict the phenomenology of 0.) 

Given this aim, we select the models defined by the action 

4 n) (p) = i|dt0 P (-<9 t 2 -^ p + i^ 

+gA LV [ dt [°°dk (JL) 4> p d t ^*(p, k) . (35) 



oo 



LV 



In have factorized out a factor of A^y so that the coupling constant g is dimensionless. 

When compared with the action of eq. ( l20l) . we have replaced the discrete index i by the 
integral over the dimensionless variable k. As recalled in Appendix A, the spectrum of the 
environment must be continuous to have proper dissipation, see discussion after eq. ( 1761) . The 
variable k can be viewed as a momentum, in the units of A^y, in a flat extra spatial dimension. 
The relationship with some Brane World Scenarios in then clear [HI [HI]. In the 'atomic' version 
of this model [El [T7] which has inspired us, the radiation field \l/ is a massless 2 dimensional 
field propagating in the dimension associated with k. 

We have also introduced an additional time derivative acting \l/ in S^^. This choice leads 
to the above mentioned 5(t — t'). Indeed, on the one hand, taking account this extra derivative, 
the continuous character of k, and the fact that g is independent of k, eq. fl2"5"]) becomes 



[d 2 t + uj 2 )(f> p = g n d t J »°(p, k, t) - g n d t J dt 1 J dkR°{t, t'; k,p) d t , (<?„0 p (t')) , (36) 

where g n = gA^y ip / ^Lv) n+l ■ On the other hand, for each 3-momentum p, \I/ = J dk^(k) is a 
massless 2-dimensional free field. In Fourier components, its retarded Green function is given 
by R°(oj, k) = l/(-(w + it) 2 + (nA LV k) 2 ). Hence R°(t,f) obeys 



X' 



d t R°(t,t') = d t I ^ rdkR°(u,k)e-^-^ = ^-A (37) 
27i A LV 



which is the required property to simplify eq. ( 1361) . 

When g is constant, the retarded Green function of associated to eq. ( 1361 ) obeys the 
following local equation 



9l 



[d 2 t + f 2 -^ + wj] G r (t, t>, p) = 6(t - , (3t 
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To make contact with Appendix A and B, let us rewrite this equation in Fourier transform: 



r 2 

[-u -ig 



A LV Alv 



(39) 



We thus see that ReS r = and that ImS r is (exactly) given by g 2 times that of eq. (|85p . 
Thus, even though we have chosen a simple form for R°(i, t'), the above action delivers the n 
dissipative behaviors of Appendix B by choosing the appropriate power of p/ Alv in S^. 

When g and u 2 are arbitrary time-dependent functions, the Fourier analysis looses its power. 
However, in time dependent settings, our G r still obeys a local equation: 



where the n-th decay rate 7 n (t) = g 2 {t)^ n is now a definite time dependent function. 
2.6 Covariant description 

We now provide the covariantized expression of the action St = J <i 3 pSV(p), where St(p) is 
given in eq. (1351) . This expression will then be used to define our theory in curved backgrounds^] 
Two steps should be done. We need to go from p considerations to a local description, and 
express the various actions in terms of the unit vector field and the spatial metric U iV . Both 
are straightforward and, in arbitrary coordinates, the action reads 



where the symbol A is the Laplacian on the three surfaces orthogonal to P. 

We have slightly generalized the action of eq. (1351) by subtracting c% U*" v to the kinetic term 
of where c|, <C 1. With this new term, the \l/(fc) are now massive fields which propagate 
with a velocity whose square is bounded by c%. In addition they now possess a well defined 
energy-momentum tensor which can be obtained by varying their action with respect to g^ v . To 
obtain the simplified expressions we have used (and shall still use), the (regular) limit — > 
should be taken. 

2 This procedure perhaps requires further explanation since ^ is not a fundamental field but "nothing more 
than a convenient parameterization of some environment degrees of freedom." First, ^ has not been introduced to 
parameterize the dissipative effects arising from any theory, but only those from theories obeying the Equivalence 
Principle, see [9] for a prototype. Second, from the point of view of cond-mat physics, it could a priori seem 
inappropriate to proceed to a covariantization, since in most situations (heat bath) there is a preferred frame 
which is globally defined. However, when the system is non-homogeneous (e.g. a fluid characterized by a non- 
homogeneous flow), low energy fluctuations effectively live in a curved geometry [28|. Moreover, in this case, 
short distance effects, i.e. dispersive (or dissipative) effects, are covariantly described when using this metric 
because they arise locally [31 [TT] . Hence the covariantization is both a necessary step to implement the EP in 
our settings and a property that emerges in cond-mat physics. 




(40) 




LV 



(41) 
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In this limit, the (free) retarded Green function of the ^ field obeys a particularly simple 
equation when expressed in space-time coordinates: 



d 



d 1 S 4 (x ti — ?/ 

dkR°(x,y;k)=l»—R°(x,y) 



(42) 



On the r.h.s, one finds the delta function with respect to the invariant measure d A x^j —g. This 
equation is nothing by the covariantized and "localized" version of eq. (j37p . Its physical 
meaning is clear. In our model, the back-reaction of <p(x) onto itself through is local. 

In spite of this, the equations of motions do not have a particularly simple form. Using the 
condensed notion = J dkty, one gets 



1 



d^g^d v ct ) {x) = gK 



1 



LV 



-(—) 



n + l 
2 



where the interacting *&(x) field is 

*(x) = *°{x) - gA LV [ d 4 y^gR°(x,y) 



/ A Tl + l \ 



-9 -LV 



(43) 



(44) 



When inserting eq. (jHJ) in eq. f|4"3]) . using eq. fj42l . one verifies that the dissipative term is local 
and first order in Fd^. Therefore, as expected, dissipation occurs along the preferred direction 
specified by the vector field Z M . 



2.7 Dissipative effects in curved background geometries 

To define a dissipative QFT in an arbitrary curved geometry, one needs some principles. From 
a physical point of view, we adopt the Equivalence Principle, or better what can be considered 
as its generalization in the presence of the unit time-like vector field Z M . We are in fact dealing 
with two (set of) dynamical fields, the field we probe, and the ^ field we do not; but also 
with two background fields g^ v and F. The Generalized Equivalence Principle means that the 
action densities of the dynamical fields be given by scalar functions (under general coordinate 
transformations) which coincide to those one had in Minkowski space time for a homogeneous 
and static F field, i.e. those of eq. fT4TT) . 

However, the densities are not completely fixed by the GEP. In this we recover what was 
obtained with the EP: For a scalar field in a curved geometry, there was always the possibility of 
considering a non-minimal coupling to gravity by adding to the Lagrangian a term proportional 
to R(p 2 . In the present case, the ambiguity is larger because P defines new scalars, the first of 
which is the expansion = V^P, where V M is the co variant derivative with respect to g^. The 
ambiguity can only be resolved by adopting some additional principle, such as the principle of 
minimal couplings which forbids adding densities containing these scalars. 

Rather than adopting it, we shall choose the non-minimal coupling so as to keep eq. ()42l) . 
i.e. so that the locality of the back-reaction effects of <p through ^ be preserved. This choice 
maintains the simplicity of the equations of motion in curved backgrounds, but is by no means 
necessary. Starting from eq. (HTl) . the locality is preserved by replacing in and Sy^ 

Fd^ k -> V&h = Fd^ k + ^® = \ (^V M * fc + V M [F* k ]) . (45) 
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To simplify the forthcoming equations, we use the fact that one can always work in "pre- 
ferred" coordinate systems in which the shift P vanishes and in which the preferred time is such 
that l — 1. (We assume that the set of orbits of Z M is complete and without caustic. In this 
case, every point of the manifold is reached by one orbit.) In these coordinate systems, it is 
useful to work with rescaled fields \I/ r = (— g) 1//4 \I/ because the above equation simplifies 

T>cs h = (-g)- 1/4 /^((-s) 1/4 **) = {-g)- 1/A l"d,% , (46) 

since 6 = (— flO -1 ' 2 f^[( — g) 1 ^ 2 ^]. Notice also that there exists a subclass of background fields 
(g,l), for which one can find coordinate systems such that both the shift l l and g m vanish. In 
these comoving coordinate systems, the above equations further simplify since only the spatial 
part of the metric matters because — g = h c where h = det(J_y). [f| 

Having chosen this non-minimal coupling, one verifies that the kinetic term of the rescaled 
fields \l/ r is insensitive to the "curvature" of both g^ u and P (when the limit c| — > is taken). 
Moreover the differential operator which acts on the retarded Green function of \l/ r in the 
equation of motion of 0, see eqs. ()4"3| 1441). is also "flat" thereby guaranteeing that the modified 
version eq. (jl2"j) still applies, that is 



T>iR°(x,y) = {-g{x)y^l»^(n° r {x,y)){-g{y)y 
1 8\x»-y») 



A 



LV 



-1/4 

(47) 



where H°(x,y) is the retarded Green function of the rescaled field \l/ r . It obeys (in pre- 
ferred coordinate systems) d t H° = <5 4 /Alv, and "defines" the retarded Green function R° = 
(— g) 1 ' A 'R° (— g) 1 ^ which is a bi-scalar. Hence the equations of motion in an arbitrary back- 
ground " tensor- vector metric" specified by the couple (g^uJ 11 ), are given by eqs. fj4"3j l4*4l with 
the substitution of eq. ( 14*51) . 

Several remarks should be made. First, from the simplified equation c^R° = 5 a /Alv it 
might seem that the background tensor metric g^ plays no role. This is not true, since it is 
used to normalize the field at every point. 

Second, in the limit c| — > 0, the (rescaled) fields define a new kind of field. They 
propagate in an effective space-time given by the time development of the 3-dimensional set 
of orbits of the Z M field. Indeed, at fixed k, ^(x) can be decomposed in non-interacting local 
field-oscillators, each of them evolving separately along its orbit. This situation is similar to 
the long wave length (gradient-free) expression of [32]. In the absence of / M , the geometry must 
be (nearly) homogenous for the action to posses this decomposition. However, when is given, 



3 It is clear that this is the case when coincides with the cosmological frame and when one uses comoving 
coordinates ds 2 = —dt 2 + a 2 dx 2 since — (1, 0). However, in certain cases one should search for the "comoving" 
coordinate system. To illustrate this point, consider the former situation in Lemaitre coordinates X = ax. In 
this case one has ds 2 = —dt 2 + (dX — Vdt) 2 , where the velocity is V = HX. The spatial sections are now 
the Euclidean space with h = 1, and the (contravariant) components of the unit vector field are I 11 = 1,V. 
To compute h c one should solve the equation of motion of comoving (free falling) observers dX — Vdt = 0, 
and use the initial position as new coordinates. This procedure is explicitely done in [12) when starting with 
Painleve-Gullstrand coordinates to describe the black hole metric and using a freely falling frame. 
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one can identify, even in non-homogeneous metrics, each space-time point in an invariant way 
by the spatial position of the corresponding "preferred" orbit at some time, and the proper time 
along the orbit (as long as Z M has no caustic). We can thus build covariant actions exploiting this 
possibility and consider fields composed of a dense set of local oscillators at rest with respect 
to Z M . The fields ^ we use belong to this class of fields. 



3 Dissipative effects in cosmology 
3.1 The action 

Our aim is to describe dissipative effects in an expanding homogeneous universe when the 
vector field is aligned along the cosmological frame, when the dissipative effects are known 
in Minkowski space, and when implementing the Equivalence Principle. 

In this case, to get the action we simply consider eq. (j4ip (with the curved metric modifi- 
cations discussed in the former subsection) in a FLRW metric. Using comoving coordinates, 

ds 2 = -dt 2 + a 2 (t)dx 2 , (48) 

the components of Z M are (1,0). To simplify the notations, we use the conformal time drj = dt/a 
and work with the rescaled fields r = a<p and \l/ r = lij 2 "^ = a?/ 2 ^! . Dropping these r indices, 
working in Fourier transform with respect to x, the non-minimally coupled action associated 
with the replacement of eq. (1451) is 



4 n) (p) = \ J dv - uftri))*. 



+ jdrjQn^^dr, Jdk^*( P ,k). (49) 

The conformal frequency u) 2 (r]) = p 2 — d 2 a/a is that of a rescaled minimal coupled massless 
field. In this expression, as everywhere in this Section, p is now the conformal (dimensionless 
and constant) wave vector. The time dependent coupling coefficient is 

g n = g a 1 ' 2 A LV (p/aA LV ) n+1 . (50) 

Its dependence in a (given that of p) follows from having implemented the Generalized Equiv- 
alence Principle (GEP) which determines the powers of a for each term in the action. As we 
shall see, this relative power guarantees that all (fi p will be damped at a fixed and common 
proper scale (in the adiabatic approximation). Had we started with an action of the type (1491) 
without relying on the GEP, the dependence in a would have been arbitrary, and the proper 
scale at which modes would have been damped would have run as well. 

In the same vein, one sees that the proper frequency of the \I/ fields stays constant. This 
follows from the GEP but also from our choice of non-minimal couplings. In this paper, we want 
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indeed to analyse the phenomenology of dissipative effects when the new degrees of freedom 
are and stay in their ground state. One could have chosen more complicated models in which 
$ is parametrically excited. This would be the case when introducing a non zero velocity 
and with minimally coupling (by dropping on the rhs of = in eq. (|45p ). However, when \1/ is 
taken massive (i.e. by restricting the k range to \k\ > 1), or discretizing k as it would be the 
case for Kaluza-Klein modes, and if Klv 3> H, the phenomenology of all these models coincide 
since the amplitude for parametric excitations will be exponentially damped. In this regime 
there is thus no gain in studying more complicated actions than that given in eq. (j49]). 



3.2 Equation of motion 



Using the results of the subsection I2.5[ the equation of motion of Heisenberg operator P is 

(d 2 . + 2 ln d v + (u 2 p ( v ) + d nln )) P = g n 3„¥°(p) , (51) 
where the decay rate in conformal time is, 

It is dimensionless, as it should be. It should be compared with the comoving frequency p to 
get the relative strength of dissipation. One obtains 

JrJjj) _ 1 / P \ 2n+l _ 1 (Pphys\2n+1 , , 

In the last equality we have re-introduced the proper momentum p p h ys = p/a. With this 
equation we verify that, at any time in an expanding universe and for every mode <p p , the 
relative strength of dissipation is simply that of Minkowski space evaluated at the corresponding 
energy scale, see eq. flHUj) . This directly results from having implemented the GEP. 

Notice however that the equation of motion in an expanding universe contains a frequency 
shift 

9n q /i„ „ \ J2 i aH \ i P 



, 2n+l 

A LV "' iy ~' c " l/ v p ' K aA LV ' 



A^9n)=P 2 (^)(-^Y n+ . (54) 



We have factorized out the unperturbed frequency square to get the relative value of the shift. 
It vanishes both when dissipation is negligible and when aH/p 1, i.e. when the expansion 
rate H is negligible with respect to the proper momentum. From this expression we can already 
conclude that it cannot play any role when the two relevant scales H and K^y are well separated, 
i.e. when 

C7 = -— «1. (55) 

Indeed when the physical momentum is high and of the order of A^y, the relative frequency 
shift proportional to a, and when the physical momentum is of the order of H (at horizon exit, 
see Figure 1), it is proportional to cr 2n+1 . 

In quantum settings however, it is not sufficient that the equation of motion possesses its 
Minkowskian form because the state of the system, and therefore the observables, might be 
affected by the combined effect of dissipative effects and the expansion rate. 
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3.3 Power Spectrum 



Let us consider the power spectrum of a scalar massless test field (which is the relevant case 
for gravitational waves and density perturbations) both in the standard free field settings and 
in the presence of dissipative effects. 



3.3.1 Free settings 

The equation of motion is the same as in eq. ( TTU]) with 

2 

n 



u 2 p =p 2 --. (56) 



For simplicity we consider the inflationary background to be de Sitter where a(rj) = —l/Hrj 
and r) < 0. At the onset of inflation, when p p h ys /H = p/a in H = p(—r] in ) — ► oo, the positive 
frequency modes are 

1 A i 



(2p)V 2 \ pr] J 

When p\rj\ = 1, the physical wave length A = a/p becomes larger than the Hubble radius. Near 
that "horizon-exit", u) 2 p flips sign, the mode stops oscillating and starts to grow like a. At late 
time with respect to horizon exit, when \pq\ <C 1, this implies that the power spectrum becomes 
a constant, as we now recall. 

When inflation lasts long enough, i.e. when the number of extra e-folds obeys N extra ^> 1, 
see eq. (pQ), all observable modes P are in their ground state at the onset of inflation (simply 
because this is the only state compatible with inflation [27]). In this case, the anti-commutator 
of the free field, see eq. ( 1321) . when evaluated at equal time r\ is simply given by 



Gf -(,, rt = W |. = _(l + _j. (58) 
Then the power spectrum of the physical (un-rescaled) field given by 

^4^=# ! < i+ « ! ). (») 

becomes constant after horizon exit. We have added a p subscript to H because in slow 
roll inflation, the relevant value of H for the p-mode is that evaluated at horizon exit, i.e. 
H p = H(tp), where t p is given by p/a = H. The above equation shows that P p acquires 
some scale dependence only through H p . Similarly the deviations from this standard behavior 
stemming from some UV modification of the theory will also depend on p through H p (and its 
derivatives). For explicit examples, we refer to [30] where the modifications of the spectrum 
stem from the fact that pr]i n is taken large but finite, and also to [3T] wherein the UV cutoff is 
endowed with a finite width. This last case bears many similarities with the dissipative settings 
we now study. 
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3.3.2 Dissipative settings 

In the presence of dissipation, the expression for G a radically differs from the above. 

At the level of the Heisenberg operator, when dissipative effects grow with the energy (as 
we suppose it is the case), the decaying solution of eq. floTj) is completely erased (unless one 
fine-tunes N extra , see eq. (pP), so as to keep a residual amplitude). That is, in inflation the mode 
operator is entirely given by its driven term, the second term in eq. ( 1261 . Then the power 
spectrum is also purely driven and given by the second term of eq. ( 1331) : 



G driuen^ p) = ^ J dr] 2 G r (r], rji, p) G r (r),T]2,p) N(r]i,T]2,p) , (60) 

where G r is the retarded Green function, solution of eq. (loTl) with 8{j] — rji) as a source, and 
where the kernel N is the anti-commutator of gnl^d^p, the source of eq. (IBTl) . 

Eq. (1601 tells us that only the state of the environment matters. In other words, because of 
the strong dissipation at early times, the power spectrum is independent of the initial state of 
<p, what ever it was. In spite of this, when eq. fl55|) is satisfied, and when the environment is in 
its ground state, the predictions are unchanged, i.e. the power spectrum obtained from G^ riven 
coincides with that obtained with G{ ree . To show this let us study G r and N. 

3.3.3 Dissipative Green functions and noise kernel 

The retarded Green function is of the form 

G r ( V , Vo ,p) = 6( V - m ) e -^>' 7(?/) x 2Im($;( V ) P M) , (61) 

where (in the under-damped regime) the modes <p p are unit Wronskian positive frequency 
solutions of eq. (lo^j) with a frequency square given by 

K") 2 = ^-7, 2 ,=P 2 d-^-|). (62) 

In this, we obtain a time-dependent version of the stationary case, see eq. flHTj) . With the 
second expression, we verify that when eq. flo^j) is satisfied, the two corrections terms are never 
simultaneously relevant, which guarantees, as we shall discuss below, that dissipative effects 
can be studied in the quasi-stationary approximation. Notice also that the frequency shift 
present is eq. (I51I) drops out from uj e " '. 
Eq. (IdTI) implies that 

G r (v,Vo,p) -> when r] -> r£ , (63) 

where r]p is the "A^y-exit" time of the p-mode defined by p/a(i]p) = A^y, or 7/p = 1/2, see 
eq. ( |53|) and the Figure. 
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Figure caption. We have represented in a log-log plot and by a dashed line the evolution of 
dn = Rh/o, the comoving Hubble radius, as a function of a, both during inflation (dn oc 1/a) 
and during the radiation era oc a). We have represented by a thick line the trajectory 
of the cutoff length d^ = of a fixed proper scale which obeys <C Rh during 

inflation. The dotted line corresponds to an intermediate fixed proper length X which obeys 
1/A LV C A < Rfl '. The vertical line represents a fixed comoving scale d p = 1/p. Below 
the cutoff length, all modes are over-damped. When a mode exits the cutoff length, it becomes 
under-damped and starts propagating. When it crosses the intermediate length X, it behaves 
as a free mode, and it gets amplified only when exiting the Hubble radius. Adiabaticity, which 
is guaranteed by l/A^v ^ Rh, guarantees in turn that, near X, modes are all born in the 
Bunch-Davies vacuum when the environment ^ is in its ground state. 

Eq. ( J63l) guarantees that no quantum coherence is left between late time 77, e.g. horizon- 
exit, and what happened for times earlier than A^y-exit. The physical implication of this is 
that, what ever transitions happened, what ever was the quantum state, no record is kept at 
late time@ Having understood the structure of the retarded Green function, let us now briefly 
consider the commutator G c (t],t]o). It is given by 



It possesses two noteworthy properties. First, because of the absolute value in the damping 
exponential term, it obeys two different local differential equations depending of the sign of 
77 — 770- When 77 > r/o, G c is an homogeneous solution of eq. (1511) . whereas when 77 < 770, it obeys 
the equation wherein the sign of the dissipative term has been changed. These two differential 

4 This behavior has been described, drawn, and sometimes named anti- dissipation, see [5J[TTJ[7J[33 ( 34, 351114). 
In fact when tracing backwards in time vacuum configurations from the freely propagating regime down to the 
dissipative one, dissipation arises towards the past. When viewed as a function of 770, G r (y, rjo) indeed increases 
to the future, thereby showing anti-dissipation. From eq. ()61j) we understand the origin of this quantum 
property (without classical counterpart). As a function of 770, first, G r obeys eq. (IST]) with the "wrong" sign for 
the dissipative term, and second, the boundary condition is a final one and not an initial one, thereby explaining 
why when 770 — > ??, the value of the current is "adjusted" to become exactly one. 




(64) 
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equations can be grouped into a single non-local equation. Second, because of this "quasi" -local 
properties, G c is independent of r]i n , the moment when was put in contact with ^f. 

These properties of G c follow from the fact that its source term, the D kernel of eq. (129j) , 
is "ultra" local for the * field of eq. (P|): 

0(771,772) = [gnd v Mm),9nd v Mri2))=9n(vi) l M^UZJ^ gn(r]2 y (65) 

This should be contrasted with the kernel N, the anti-commutator, which is non local. Indeed, 
in the vacuum, it is given 

N(?7l, 7/2) = 9n(Vl) 7° ( ?i )0( ?^ 9n(V2) , (66) 

iTA LV (ti - t 2 y 

where l/(ti — t 2 ) 2 should be understood as the derivative of the principal part of l/(ti — £2)- 
One gets a simple expression in terms of the proper time t because the proper frequency of 
is constant in our model (149 p . We also remind the reader that it is only in the high temperature 
limit that iV would become proportional to 5(ti — £ 2 )- 

In this paper we shall compute the power spectrum when the are all in their ground state. 
The reason for this choice is as follows: As pointed out in the footnote EJ the fields should 
be conceived as parametrizing some degrees of freedom in a fundamental theory obeying the 
Equivalence Principle. In these theories, when inflation lasts long enough, all relevant degrees 
will be in their ground state, as is the case when dealing with free modes, see eq. We see 

no reason why dissipation could possibly invalidate this conclusion. 

When the *f?k are all in their ground state, the kernel N is non-local. Therefore the anti- 
commutator G a obeys a non-local equation: eq. ( l5Ti) with N as the source. The physical 
meaning is that the environment is driving (p p in an non-instantaneous way (this is unavoidable 
in the vacuum since only one sign of the frequency is present). The consequence is that the value 
of power spectrum depends (to a certain extend) on the history of the combined evolution of 
</> + \E r . The mathematical implication is that an exact calculation of G a is effectively impossible. 

Nevertheless when H/Alv 1, these non local effects give only but sub-leading (non- 
adiabatic) corrections, because the evolution of <p + ^ consists in a parametric (adiabatic) 
succession of stationary states ordered by the scale factor a. 

3.3.4 Scale separation and adiabatic evolution 

Given the importance of this result, we shall spend some time to explain its root and its 
validity. We proceed in two steps. We first show that adiabaticity is sufficient for obtaining 
the standard power spectrum, and then show that scale separation, eq. ( 1551) . is sufficient to 
guarantee adiabaticity, thereby generalysing the results of [8]. 

Adiabaticity means that the evolution proceeds slow enough for inducing no (non-adiabatic) 
transition (in molecular physics, they are called Landau-Zener transitions, in the present cos- 
mological context they correspond to pair creation of <fi quanta). When this is the case, the 
observables take, at any time, the value they have in the corresponding stationary situation. 
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As of retarded Green function and the commutator, this principle does not bring any new 
information because in the WKB approximation, they are local functions since these 2pt func- 
tions obey local equations. For the anti-commutator, adiabaticity guarantees that when rj and 
rf are close (in the sense that 1 — a(r])/a(i]') <C 1), its value is well approximated by 

Gt i ^(r ] ,r]',p)^G s : atio (r ] -r ] '-,u p (a),g p (a)) = J ^ e^'> G.fa u p (a), g p (a)) , (67) 

where G a (ou;ui p (a), g p (a)) is the Fourier component calculated with the values of the frequency 
ujp{a) and the coupling g p (a) both evaluated with a = a(r]). In Appendix A, these Fourier 
components have been algebraically solved for all frequencies and all couplings, see eq. (181j) . 
Moreover, since by hypothesis, we are in the vacuum, eq. ( 1791) also applies. In other words, the 
value of G a follows from that of the commutator G c . This is sufficient to guarantee that when 
the mode (j) p becomes free, i.e. much after A^y-exit but before horizon-exit (H <C p/a <C A LV ), 
eq. ( !84l) applies. Thus, irrespectively of what was the coupling with the environment, 

G d r en {^ ri,p) - Gl ree ( V , r],p) = -i— . (68) 

2u p {a) 

With this equation we reach our first conclusion: in the adiabatic approximation and when 
the environment is in its ground state, the modes <f) p are born in the usual in-vacuum (Bunch- 
Davies) as they become, one after the other, freely propagating after A^y-exit, see the Figure. 
It now behooves to us to show that scale separation guarantees that eq. flBTl) provides a valid 
approximation for Q^ rmen before horizon-exit, i.e. when the second term in the parenthesis in eq. 
( 1621 can still be neglected. To this end, we provide an upper bound for the probability amplitude 
of getting a non-adiabatic transition in the under-damped regime, i.e. H -C pj 'a < Alv- This 
amplitude is governed by the relative frequency change 



[UJ p 



eff\2 ^2)3/2 V )<2 _ 2 \ 3 /2 

p ) \y in) \rphys i pays/ 



where p p h ys = v/ a an d Iphys — l/ a are the physical (proper) momentum and decay rate. 
Therefore, going backwards in time from the free regime up to p phys > ^lph ys (i- e - Pphys < A LV ), 
the non-adiabatic parameter steadily grows but stays bounded by 

< 3(2n + 1) a (^y n+1 < 3(2n + 1) a < 1 . (70) 

This guarantees that the amplitude for the system to jump out of the ground state is bounded 
by a (up to an overall factor which plays no role). There is no need to study the stability 
of the ground state in the transitory regime from under-damped to the over-damped modes 
because what ever transitions happened their residual impact at late time would be suppressed 
by a factor e~$ ldt ~ exp(— l/a(2n + 1)) 1. This completes the proof that scale separation 
guarantees adiabaticity, in agreement with the conclusion reached in [HI EI] . 



5 This behavior is in accord with the procedure of [29l [30l [31] wherein the initial state is imposed on modes 
when their proper momentum reaches a given value. This radically differs with the "initial slice" approach [36 
wherein the state is imposed at the same time on all modes, irrespectively of the value of their momentum. 
From what we see, the Equivalence Principle allows only the first to be dynamically realized. 
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3.3.5 Perspectives: Beyond the adiabatic approximation and scale crossing 

It is first interesting to point out that dissipative models are more "robust" than dispersive 
models in that the "driven" power spectrum determined by eq. floTJ]) is well defined even 
when H > A^y, contrary to the fact that most dispersive models make sense only when scale 
separation (H <C A LV ) is realized. Therefore, with dissipative models, one can study scale 
crossing, i.e. the crossing of H p through Alv from above to below during slow roll inflation. 

The power spectrum has been studied in this regime when using a particular Brane World 
scenario in [10J. We agree with most of the conclusions but not with that reached after Eq. 
(35) according to which "towards the end of inflation, perturbations on the brane are the sum 
of two independent fields" . We do not agree because in those settings as well the field operator 
is purely driven. Hence we do not see how to decompose it as a superposition of two commuting 
operators. Moreover, with the approximation used in [ID] , it seems that the ETC, eq. ( TTTT) . is 
violated by a amount proportional to the reported change of the power spectrum. 

In any case, the calculation of the modifications of the power spectrum introduced by 
dissipation is difficult because the non-local properties of eq. flDD]) cannot be neglected. We are 
not aware of any analytical treatment, and we are currently studying the modifications using 
numerical techniques [38] . 

Added note. Since the present paper was submitted, we have obtained several results [38J. 
In particular, the leading modification of the power spectrum with respect to the standard 
result in the regime a = H/A <C 1, i.e. the signature of UV dissipation, scales with a power of 
a which is equal to that of P/A in j/p, see eq. ( 1531) . In this respect, dissipative models behave 
like dispersive models where the leading modification scales with a power of o which is that 
of the first non-linear term in the dispersion relation, i.e. 2n using the parameterization of eq. 
([2]), see [39]. In the opposite regime, when H > A, no universal behavior is obtained. 

4 Conclusions 

In this paper we have obtained the following results. 

First, we have provided a class of unitary models defined in Minkowski space which are 
characterized by the power of p/A^y which weighs the growth of dissipation in the preferred 
frame, see eq. ( [85]) . Unitarity is achieved by introducing a dense set of additional fields 
\I/ which induce dissipation through interactions with <fi. In Appendices A and B we have 
given a thorough analysis of the Green functions of these Gaussian models which cover the 
phenomenology of UV dissipative effects in Minkowski space at the level of 2pt functions. 

Second, amongst the various actions delivering the same stationary phenomenology, we have 
selected one which gives rise to a local differential equation for the retarded Green function 
of (ft, see eqs. ([361 1371) . Using the Equivalence Principle, we have extended this class of 
models to arbitrary curved backgrounds in subsection 12.71 thereby allowing to confront the 
trans-Planckian question of inflationary cosmology and black hole physics. 

Third, we have applied our dissipative models to inflationary cosmology. At early times, 
hence at very high energy, the dissipative effects are so strong that all information about the 
initial state of the is erased, see eq. ( [63]) . Nevertheless, when the UV scale Alv is much 
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larger than the Hubble parameter, we have demonstrated that the standard expression of the 
power spectrum is found when the environment is in its ground state. The reason for this is 
the following: even though the field oscillator <p p is purely driven by *& p , i.e. it is given by the 
second term eq. fl2"B"j) . as its proper momentum p/a redshifts under A LV , the composite operator 
behaves as if it were a free operator, see eqs. fl3"Tj 1831) . thereby guaranteeing eq. (jS8|) . 

Let us also note that our models can be used for phenomenological purposes in the sense 
of pQ, and that we are planning to apply them to study Hawking radiation in the presence of 
dissipation. We are presently completing the calculation of the power spectrum beyond the 
adiabatic approximation so as to determine the signatures of dissipation [38] and to compare 
them to those of dispersion [3_H]. Finally it would be interesting to compute the VEV of the 
stress-tensor of the fields in non-trivial backgrounds. This would allow to take into account 
the back-reaction on the cosmological metric engendered by the fields St and / M [4*01 I4T1 132] . 

Acknowledgments. 

I am grateful to Dani Arteaga and Enric Verdaguer for common work allowing me to deepen 
my understanding of dissipative effects. I am also grateful to Julian Adamek, David Campo, 
Ted Jacobson, Jean Macher, Jens Niemeyer, and Valeri Rubakov for interesting discussions. I 
wish to thank the organizers of the workshops on " Micro and Macro structure of spacetime" 
held in Peyresq in June 2005, 2006, and 2007 where this work has been presented and discussed. 
This work has been supported by the Agence Nationale de la Recherche (projet 05-1-41810). 



5 Appendix A : 

Stationary states and vacuum 2pt functions 

In this Appendix, we recall the (well-known) relationships between G c , G a which always hold 
in stationary states. In these states, G c , G a and the kernels D, N are related to each other in 
a universal way, generally referred as a Fluctuation-Dissipation relation. We explain its origin 
and its physical implications in the present context. We start with the most basic object: the 
retarded Green function G r . 

5.1 The retarded Green function 

The Fourier transform of eq. (I25p gives 

(-u 2 + ^)0 P H = E igi (p)*?(p,u;) + E^ 2 (p)i?°(u;;p)0 P H , (71) 

where 

i2?( W ;p)=(-(o; + ie) a + ^(p)))- 1 ) (72) 

is the Fourier transform (defined as in eq. (1171) ) of the retarded Green function of As usual, 
its retarded character is enforced by the imaginary prescription of the two poles to lay in the 
lower half plane (e > 0). The solution of the above equation is 

<P p (tu) = cf> d p {u) + G r {u,p)Z igi (p)*°(p,u) , (73) 
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where the Fourier transform of the retarded Green function of <j), the solution of eq. ( J271) . always 
takes the form 

G r (u,p) = (-(u + ief + p 2 + Z r (uj,p)y 1 . (74) 

All effects of the coupling to the \I/j's are thus encoded in the (retarded) self-energy T, r (u,p). 
For Gaussian theories, it is algebraically given by 

Z r (uj,p) = -i; i g 2 (p)R°(u J ,p). (75) 

The dissipative effects are governed by the imaginary part of E r (a;,p). In the present case, one 
has 

2ImZ r (u,p) = -Z ig 2{p)G° Cti {u,p) = —D(uj,p) . (76) 

To get the first equality we have used the fact that in stationary states the retarded Green 
function and the commutator are related by 2lmG r (uj) = G c (oj) for all degrees of freedom, free 
or interacting. In the second equality, we have introduced D(lv), the Fourier transform of the 
kernel of eq. (IHUj) . 

Several observations should be made here. First, from eq. (1721 . we obtain that D{uj) is 
proportional to Y>ig 2 5(u) — Qi). Therefore there is no dissipation for lower frequencies than the 
lowest value of This simply follows from energy conservation. Second, to obtain "true" 
dissipation, D(u,p) should be a continuous function and not a sum of delta. This can only 
happen when the form a dense ensemble. In Section 12.5} we shall thus replace the discrete 
sum on i by an integral over a continuous variable, k. We shall not consider the discrete cases 
even though these could display interesting properties. Third, from a phenomenological point 
of view, only D(u,p) matters. Hence we cannot disentangle the spectrum of the environment, 
which is given by R°(u,p), from the coupling strength gfip). This is a good thing, because 
when working in time-dependent settings, we shall exploit this equivalence to chose the simplest 
model of \l/j's which gives the kernel D(u),p) we want. 

It is also worth noticing that the dispersive (real) effects are not directly related to D (or 
N). These are governed by the even part of H r (u,p) which is given by 

ReE r (u,p)= r , (77) 

where the integral should be understood as a principal value. This integral relationship explains 
why one often founds that dispersive effects appear before dissipative effects (for increasing u). 
We also learn that the dispersive models studied in the literature violate the above relation 
since they assume ReS r ^ and ImS r = 0. Therefore these models are incoherent and cannot 
result from dynamical processes. 

5.2 Fluctuation-Dissipation relations and vacuum self-energy 

In this subsection, we derive the relationships between G c , G a and T,p which exist in the true 
(interacting) ground state. 
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In interacting theories, the only stationary states are thermal states, see e.g. [37] . In these 
states, the Fourier transform of D and N are related by 

N(u) = D(cu) coth(/3w/2) , 

= D(u) sign(u;) [2n(\w\) + 1] . (78) 

In the second line, n(u) is the Planck distribution. It gives the mean occupation number of 
quanta as a function of the frequency (measured in the rest frame of the bath). The above 
relation directly follows from the fact that the individual commutators and anti-commutators 
of the free fields obey this relation, as any free oscillator would do. It implies that the 
Fourier transform of G c and G a are also related by 

G a {u) = G c {u) sign(w) [2n{u) + 1] . (79) 

It should be stressed that this equation is exact, i.e. non-perturbative, and valid for all theories, 
Gaussian or not. (It indeed directly follows from the cyclic properties of the trace defining 

G^(M0 =Tr[e-^(tW)]). 

For Gaussian models, there exists an alternative proof of eq. ( 1791) . It suffices to note that 
in steady states the decaying terms of eqs. ( 1291) and (1331) play no role, and that the Fourier 
transform of the driven terms are respectively given by 

G c {u) = \G r (u)\ 2 D{u) , (80) 
G a (u) = \G r (u)\ 2 N(u), (81) 

since the Fourier transform of the retarded Green function obeys G r (u) = G*(—u), see eq. 
( 1731) . Irrespectively of the complexity of G r , i.e. irrespectively of the functions gi(p), £li(p) and 
the set of the ^ fields, G c and G a are thus related to each other by the FD relation (179"]) . 

These universal relations will be relevant for inflationary models wherein only the ground 
state contributes. In particular, they imply that in the true vacuum, i.e. when n(u) = 0, G c 
and G a are exactly related by G a (u) = G c {uj) sign(u;). Hence the Wightman function 

G w = l -(G c + G a )=G c 6(uj) ) (82) 

is determined by the commutator and contains only positive frequency, as in the free vacuum. 
Equations ( IHUllHTi) also allow to compute the vacuum self-energy of the Feynman Green function. 
For Gaussian models it is given by 

2ImE^(a;,p) = —D(u,p) sign(cj) . (83) 

With the last equality we recover the fact that in the vacuum, it is sufficient to consider Feynman 
Green functions. In non-vacuum states, and in non-stationary situations, this is no longer true, 
thereby justifying the use of the in — in machinery [20] (the Schwinger-Keldish formalism). 

Before specializing to a specific class of models giving rise to dissipation at high frequency, 
we make a pause by asking the following important question: What should be known about 
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the fields to get eqs. ( J80l IHTT l82l) ? We have proven that it is sufficient for the \tVs to be 
canonical fields, but is it necessary ? 

The answer is two fold. On one hand, the cannot be stochastically fluctuating quantities 
(i.e. commuting variables) because this would lead to a violation of eq. f!75|) that would imply 
the violation of eq. (JT^D and the ETC eq. (IHTT)H They cannot be either a combination of 
quantum and stochastic quantities because this would still lead to a violation of the ETC. 
Hence they must be built only from quantum (canonical) degrees of freedom. 

On the other hand, the \!Vs can be composite operators, i.e. polynomials of some (unknown) 
canonical fields. Indeed, their commutators would still be all related to their anti-commutators 
by the FD relation eq. (ITS]) , and this even though they depend non-linearly on n(u) in non- 
vacuum states. The difference with Gaussian models is that these non-linear operators posses 
non- vanishing higher order correlation functions. Hence, the self-energies S r , will be series 
in powers of gf, and not just a single quadratic term as in eq. ( 1751 1831) . Nevertheless these 
higher loops corrections preserve the validity of eq. (|82|) in the ground state, as well as that 
of eqs. (ISUl IHTT) in any thermal state, when properly understood, i.e. with D now defined by 
-2ImS r (as the effective dissipation kernel), and N related to it by the FD relation. 

In brief, we have reached/recalled the following results. Firstly, the q-number combination 
\I/ = Ej(yfj\l/i, the fluctuating source term of 0, must obey the FD relation ( 17H1) . This can 
either be postulated, or better, be viewed as resulting from the fact that ^ is entirely made 
out of quantum degrees of freedom. Secondly, to lowest order term in g, the self-energy can 
be obtained by treating $ as a Gaussian variable, what ever its composition may be. Thirdly, 
when dealing with non-Gaussian theories, once having computed £ r (u;), the resulting equations 
for the 2pt functions have the same structure and the same meaning as in Gaussian theories, 
with D replaced by -2ImS r . Therefore, the phenomenology of two-point functions is entirely 
covered by Gaussian settings. 

5.3 The double limit: g 2 T — > oo followed by g 2 — > 0. 

To perform a phenomenological analysis of dissipation, we need to understand how the theory 
behaves in transitory regime from dissipative to free propagation. Similarly, to study primordial 
spectra in inflation or Hawking radiation, we also need to understand how free motion emerges 
as the proper frequency get red-shifted. It is therefore relevant to study the behavior of the 
two-point function in the following double limit. 

One first takes g 2 T — > oo, where T = t — ti n , ti n being the moment when the interactions 
are turn on, and t the moment when one studies the field properties. This limit implies that 
the decaying term in eq. (|2"oT) plays no role. Therefore, near time t, the Heisenberg field </>(t) is 
a composite operator which only acts in the Hilbert space of 

Secondly, one considers the "free" limit g 2 — > of that composite operator. One could 
naively conclude that G c and G a of eqs. ( ISUl IHTT) would vanish since both D and N are 

6 This constitutes the simplest proof that it is inconsistent to couple quantum variables to stochastic (or 
classical) ones. If one does so, the ETC of the dressed quantum variables will always be dissipated after a time 
of the order of 7 -1 . One can therefore view the experimental evidences for the ETC of some degrees of freedom 
as a very strong indication that all dynamical variables in our world are quantum mechanical in nature. 
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proportional to g 2 . However, this is not the case, because the common prefactor, |G>| 2 , is 
singular in this second limit. In fact, one verifies that it scales in 1/g 2 in such a way that, in 
the (interacting) vacuum, one always recovers 



This is the standard vacuum fluctuations of a free massless mode of momentum p. 

Two important lessons have been reached. First we learned is that even though acts only 
on the ^-Hilbert space, when g 2 — > 0, it behaves as if it were a free mode possessing its own 
Hilbert space, with no reference to ^-dynamics. Secondly, the quantum state in this would be 
Hilbert space is still exactly that of Therefore, in stationary situations, the only "souvenir" 
kept by the composite operator is the equilibrium distribution n(u) inherited from its parents. 

Let us now emphasize that the above limit is relevant for non-Gaussian models as well. 
Indeed, in the limit g 2 — > 0, there will always be a value of g 2 sufficiently small that the 
model can be well approximated by a Gaussian model. Therefore the behavior of the 2-point 
functions in the transitory regime from dissipation to free propagation can be analyzed by 
studying Gaussian models (at least in the quasi-static limit). 



6 Appendix B: 

Dissipative effects above A^y The Phenomenology 



We now have all the tools to construct models giving rise to dissipation in the vacuum above a 
critical energy scale A^y. In this Appendix we work from a purely phenomenological point of 
view, and provide the class of dissipative models wherein the imaginary part of the self-energy 
is governed by a single term, in analogy with the dispersion relations of eq. (El). 

From a phenomenological point of view, if one considers only stationary situations (i.e. static 
metrics and stationary states), one can simply choose the function D{oj,p) entering eq. ( I76j ) 
and eq. (1HU1) as one wishes. There is indeed no restriction on D(u,p) besides its constitutive 
properties, namely being odd in u and giving rise to poles in G r all localized in the lower half to 
plane. In this we have reached our first aim, namely identify how to generalize the free settings 
so as to incorporate some arbitrary dissipative effects. 

We can thus consider the dispersive models which correspond to those defined by eq. (j2J). 
They are characterized by a single term giving rise to dissipation above A^y. In the vacuum, 
they are fully specified by the imaginary part of the (retarded) self-energy 







In these models, the decay rate (inverse life time) on the mass shell is 
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To verify it, assuming that Re£ r = 0, the two poles of G r (u) in eq. (1741) are located in 



w ± (p) = ± v /^T^_, 7 . ( 87 ) 

From this, by inverse Fourier transform G r {uj), one obtains that the decay rate is indeed 7 in 
the under-damped regime, for > 7 2 . In the overdamped regime, for 7 2 > cu 2 , the decay rates 
of the two independent solutions of G> = are r± = 7 ± y/7 2 — cu 2 . 

One thus have the following behavior as p grows. For p <C Alv, u± — P, and one has 
a free propagation which is slightly damped with a life time in the units of the frequency 
given by (A LV /u) n+1 3> 1. In the opposite regime of high momenta p ^> A LV , deep in the 
overdamped regime, the two roots u± are real and the notion of propagation (in space-time) 
is absent. In anticipation to what will occur in inflation (or in black hole physics), we invite 
the reader to study the migration of the poles of G r when extrapolating backwards in time a 
mode, i.e. as p increases. (Remember that the physical momentum of a mode in cosmology is 
Pphys{t) oc p/a(t) where p is the norm of the conserved comoving wave vector (near a black hole 
one finds p(r) oc u/x where x = r — r$ is the proper distance from the horizon measured in a 
freely falling frame, and u the conserved Killing frequency.)) 

One could of course generalize the above class by considering in eq. ( 1851) polynomials in 
p dimensionalized by different UV scales. However, unless fine tuning, the phenomenology of 
the transition from the IR dissipation-free sector to the dissipative sector will be dominated a 
single term. One should also consider the possibility that ImS strictly vanishes below a certain 
frequency Qi, as this would be the case when ever the spectrum of the \I> fields possesses this 
gap, see the remarks after eq. f l76l ). 

Having the phenomenology of dissipative and unitary models under control (with dispersive 
and dissipative effects related by Kramers relation, eq. ([77]) ) one could confront particle and 
astro-particle physics data and put lower bounds on A LV for each n, in analogy with what was 
done for (pure) dispersion in pp. 
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